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Abstract 


Sulfur  mustard  (HD)  is  an  alkylating  agent  that  reacts  rapidly  with 
macromolecular  targets  resulting  in  the  formation  of  stable 
adducts  providing  depots  for  markers  of  exposure.  The  purpose  of 
this  study  was  to  validate  an  analytical  procedure  for  detection  of 
HD-plasma  protein  adducts  and  to  establish  the  utility  of  the 
method  in  an  HD  rat  inhalation  study.  Calibration  curves  were 
prepared  in  human  and  rat  plasma  at  six  levels  of  HD  (12.5  to 
400nM).  Correlation  coefficients  for  the  mean  data  were  0.9987 
for  human  and  0.9992  for  rat  plasma.  The  percent  coefficient  of 
variation  (%CV)  derived  from  the  mean  concentration  data  ranged 
from  0.53  to  14.1  %  in  human  (n  =  5)  and  0.57  to  10.63%  in  rat  (n 
=  6)  plasma.  Intraday  and  interday  precision  and  accuracy  studies 
were  conducted  at  three  concentration  levels  (25, 1 50,  300nM)  to 
represent  low,  medium,  and  high  concentrations  of  HD  relative  to 
those  employed  in  the  calibration  curve.  Precision  and  accuracy 
were  assessed  by  determining  %CV  and  %  error,  respectively.  For 
intra-  and  interday  studies,  the  %CVs  and  absolute  %  errors  were 
less  than  15%.  The  limits  of  quantitation  were  20.88nM  for  human 
and  16.73nM  for  rat  plasma.  In  animal  studies,  rats  received 
nebulized  HD  at  six  doses.  The  data  indicate  a  dose-dependent 
relationship  between  maximal  plasma  concentrations  and  dose 
administered  (Rl  =  0.9728).  Results  from  this  study  indicate  an 
accurate,  precise,  and  sensitive  method.  The  method  was  useful  in 
determining  plasma  protein  adduct  formation  in  a  rat  inhalation 
model. 


Introduction 

Exposure  to  vesicating  agents  such  as  2,2'-dichlorodiethyl 
sulfide  (sulfur  mustard;  HD)  is  a  concern  of  both  civilian  and 
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military  populations.  Sulfur  mustard  is  an  alkylating  agent 
that  reacts  rapidly  with  water,  peptides,  proteins,  ribonucleic 
acid,  and  deoxyribonucleic  acid  (1).  Reactions  with  macro- 
molecular  targets  can  result  in  the  formation  of  stable  adducts, 
which  provide  depots  for  markers  of  exposure  (2-4).  Theoret¬ 
ically,  macromolecular  markers  may  he  present  until  elimi 
nated  by  metabolic  processes.  These  sources  can  potentially  be 
leveraged  as  methods  to  verify  exposure  to  the  agent  at  rela¬ 
tively  long  periods  following  the  event  (2-4).  Reviews  on  assays 
designed  to  detect  exposure  to  chemical  warfare  agents  in¬ 
volving  free  metabolites  as  well  as  macromolecular  adducts 
have  been  published  (2-5).  For  the  latter  group,  assays  in 
elude  the  immunoassay  of  HD  bound  to  DNA  (6),  the  gas  chro¬ 
matography-mass  spectrometry  (GC-MS)  assay  for  the 
N-terminal  valine  adduct  of  hemoglobin  (7-10),  a  GC-MS 
assay  of  thiodiglycol  (TDG)  from  skin  keratin  (9),  the  liquid 
chromatography-MS-MS  assay  of  the  cysteine  adduct  in  al¬ 
bumin  (11),  and  an  immunological  method  for  the  detection 
of  mustard  adducts  to  skin  tissue  (12). 

The  studies  in  this  current  report  describe  validation  and 
further  use  in  an  HD  inhalation  study  of  a  previously  reported 
method  for  detecting  HD  exposure  (13).  The  technique  lever¬ 
ages  interaction  of  HD  with  plasma  proteins  to  form  stable 
adducts.  The  analytical  methodology  involves  isolation  of 
plasma  proteins,  base  catalyzed  cleavage  of  the  adduct  fol¬ 
lowed  by  chemical  derivatization  and  GC-MS  assay  of  the 
derivative.  Similar  to  skin  keratin,  plasma  proteins  contain 
numerous  free  carboxylic  acid  groups  from  aspartic  acid  and 
glutamic  acid  that  can  be  alkylated  by  the  electrophilic  HD  to 
give  hydroxyethylthioethyl  (HETE)  esters  (Equation  1).  These 
interactions  provide  a  relatively  high  quantity  of  plasma  pro- 
tein-HD  adducts  when  compared  with  adducts  formed  at  the 
four  terminal  valine  sites  in  globin  or  the  single  cysteine  site 
in  albumin  (2). 

The  interaction  of  HD  with  proteins  and  the  subsequent 
steps  in  the  analytical  method  are  based  on  the  key  reactions 
shown  in  Scheme  1.  Formation  of  the  protein-HETE  ester 
that  occurs  in  vivo  is  described  in  Equation  1.  These  esters  can 
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Equation  1 

Protein-COO  +  S(CH2CH2C1)2  +  H20  -»  Protein-COOCH2C'H2SCH2CH2OH  +  HC1  +  Cl 
(protein)  (IID)  (protein-HETE  ester) 

Equation  2 

Protein-COOCH2CH2SCH2CH2OH  3  NaOH  -  S(CH2CH2OH)2  +  Protein-COONa 
(protein-HETE  ester)  (TDG) 

Equation  3 

S(CH2CH2OH)2  +  2  CoFjCOCI  -»  (C6F5COOCH2CH2)2S  +  2  HCI 
(TDG)  (TDGPFB) 

Scheme  t.  Key  reactions  of  HD  with  proteins. 


be  readily  cleaved  or  hydrolyzed  with  dilute  base  to  yield  TDG 
(Equation  2).  The  TDG  can  be  derivatized  with  pentafluo- 
robenzoyl  (PFB)  chloride  to  give  the  bis(pentafluorobenzoyl) 
ester  of  thiodiglycol  (2,2'-thiobisethanol  dipentafluoroben- 
zoate,  TDGPFB)  (Equation  3)  and  analyzed  by  a  GC-MS  neg¬ 
ative  ion  chemical  ionization  (NICI)  in  the  selected  ion 
monitoring  mode  (9,14,15). 

Although  this  method  has  been  utilized  in  animal  studies 
and  in  an  accidental  human  exposure  to  HD,  it  has  not  been  fully 
characterized  (13,16).  The  purpose  of  this  study  was  to  1.  char¬ 
acterize  and  validate  the  plasma  protein  adduct  assay  in  plasma 
and  2.  demonstrate  the  utility  of  the  method  in  a  procedure  in¬ 
volving  the  inhaled  administration  of  HD  as  a  nebulized  aerosol 
directly  into  the  lung  of  a  rat.  Resul  ts  from  this  study  indicate  an 
accurate,  precise  and  sensitive  method.  The  method  was  useful 
in  determining  plasma  protein  adduct  formation  in  a  dose-de¬ 
pendent  fashion  in  a  rat  inhalation  model.  Herein,  we  report  the 
details  of  the  method  validation  and  its  use  in  vivo. 


Experimental 

Materials 

Ethyl  acetate,  acetone,  sodium  bicarbonate,  anhydrous 
sodium  sulfate,  pyridine,  and  PFB  chloride  were  obtained  from 
Sigma-Aldrich  (St.  Louis,  MO)  and  used  without  purification. 
Sulfur  mustard  was  obtained  from  the  U.S.  Army  Edgewood 
Chemical  Biological  Center  (Aberdeen  Proving  Ground,  MD); 
purity  as  determined  by  nuclear  magnetic  resonance  spec¬ 
troscopy  was  97.5%.  The  internal  standard  octadeuteroth- 
iodiglycol  (TDG-d8)  was  obtained  from  Ash  Stevens  (Detroit, 
MI).  Bond  Elut  silica  100-mg  Si  solid-phase  extraction  (SPE) 
cartridges  (\arian,  Harbor  City,  CA)  were  conditioned  with  1 
mL  of  ethyl  acetate  before  use.  Pooled  human  and  rat  plasma 
derived  from  sodium  heparin  anticoagulated  whole  blood  was 
obtained  from  Bioreclamation  (Hicksville,  NT). 

Animals 

Animals  were  maintained  and  used  under  a  program  ac¬ 
credited  by  the  Association  for  Assessment  and  Accreditation 
of  Laboratory  Animal  Care  International  (AAALAC  Interna¬ 
tional).  The  experimental  protocol  was  approved  by  the  Animal 
Care  and  Use  Committee  at  the  United  States  Army  Medical  Re¬ 
search  Institute  of  Chemical  Defense,  and  all  procedures  were 
conducted  in  accordance  with  the  principles  stated  in  the 
Guide  for  the  Care  and  Use  of  Laboratory  Animals  (National 


Research  Council,  1996)  and  the  Animal  Wei 
fare  Act  of  1966  (P.L.  89-544),  as  amended  in 
an  AAALAC  International  accredited  facility. 
Male  Sprague-Dawley  rats  (240-270  g)  were 
used  for  all  inhalation  experiments.  Animals 
were  anesthetized  intramuscularly  (im)  with 
the  co-administration  of  xylazine  (10  mg  kg) 
and  ketamine  (90  mg/kg).  Following  anes¬ 
thesia,  rats  were  intubated  with  glass  Pasture 
pipette  endotracheal  tubes  that  were  trimmed 
to  7  cm  in  length  and  fire-polished  and  taped 
in  place  to  the  rostrum  of  the  rat  using  surgical  tape.  Rats  were 
then  placed  in  a  glovebox  exposure  chamber  (Baker,  Sanford, 
ME)  that  was  under  a  slight  negative  pressure  (-1.0  to  -0.5  cm 
H20).  The  endotracheal  tube  was  connected  to  a  nebulizer 
system  (Aeroneb.  Aerogen,  Dangan,  Galway,  Ireland),  which 
were  in  series  between  the  rat  and  the  ventilator.  Sulfur  rnus 
tard  was  delivered  by  the  ventilation  action  of  a  rodent  venti 
lator  (Harvard  model  883)  that  was  set  to  a  respiratory  rate  and 
tidal  volume  of  60-100  breaths/min  and  2.5  mL,  respectively. 
The  HD  was  then  added  to  the  nebulization  mixture,  which 
consisted  of  25%  ethanol  and  75%  saline.  The  nebulizer  con¬ 
sistently  delivered  a  particle  size  of  <  3  pm  at  a  rate  of  0.5 
mL/min  as  determined  by  an  aerodynamic  particle  analyzer 
(TSI  model  3321,  Shoreview,  MN).  The  aerosol  stream  was 
driven  into  the  rat  through  a  two-way  non-rebreathing  (Hans 
Rudolph,  Kansas  City,  MO)  valve.  Exhaled  HD  exited  the  rat 
through  the  outflow  port  of  the  Rudolph  valve,  into  a  flask  con¬ 
taining  0.5%  bleach,  and  then  through  a  charcoal  filter  prior 
to  HEPA  filtration.  The  entire  exposure  circuit  was  a  closed 
system. 

To  verify  III)  vapor  concentrations,  samples  were  obtained 
using  a  100-pL  gas-tight  syringe  at  the  entrance  of  the  glass 
endotracheal  tube  approximately  12-15  cm  from  the  outflow 
of  the  nebulizer.  The  concentrations  were  determined  using 
Minicams IM  (O.I.  .Analytical,  CMS  Field  Products  Group,  Birm 
ingham,  AL)  chemical  agent  monitors.  The  sample  was  loaded 
into  a  10-ft  heated  sample  line  that  ran  through  the  glovebox 
wall  and  was  pulled  directly  onto  a  preconcentrator  tube  (100- 
mm  length,  3-mm  i.d.,  15-mm  bed  of  Tenax)  within  the  Mini¬ 
cams.  Separation  was  achieved  with  a  DB-1  column  (15-m, 
0.32-mm  i.d.,  5-pm  film  thickness,  J&W  Scientific,  Folsom,’ 
CA)  with  a  helium  flow  rate  of  40  mL/min.  Detection  was  ac¬ 
complished  with  a  flame  photometric  detector.  Samples  were 
quantified  against  HD  calibration  curves  generated  from  stan¬ 
dards.  The  programmed  cycle  time  between  sample  collection 
and  analysis  was  7  min,  which  allowed  the  analysis  of  two 
samples  within  the  10-min  exposure  time  frame. 

For  this  dose-response  study,  rats  received  a  single  exposure 
to  six  levels  of  HD  aerosol  at  the  following  doses  25,  75, 175, 
225,  275,  and  300  pg  HD/min  for  10  min.  The  total  dose  de¬ 
livered  to  the  animals  over  10  min  was  250,  750, 1750,  2250, 
2750,  and  3000  pg  HD.  Animals  were  deeply  anesthetized  with 
the  im  co-administration  of  xylazine  (10  mg/kg)  and  ketamine 
(90  mgkg)  and  then  euthanized  via  exsanguination  at  0,  0.5, 

1, 3, 6,  and  24  h  after  exposure  to  obtain  whole  blood  samples 
(4  6  animals  per  time  point).  Protein  obtained  from  plasma 
( 1 .0  mL)  was  prepared  as  described  in  Sample  preparation. 
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Sample  preparation 

To  obtain  plasma,  whole  blood  was  centrifuged  at  2000  x  g 
for  30  min  at  5°C;  the  plasma  supernatant  was  then  removed 
for  protein  precipitation.  Protein  from  plasma  samples  ob¬ 
tained  from  exposed  animals  or  spiked  standards  was  isolated 
and  prepared  as  described  by  Capacio  et  al.  (13)  with  modifi 
cations.  Ten  milliliters  of  acetone  was  added  to  1  m.L  of  plasma, 
which  was  then  placed  on  a  shaker  for  10  min.  Following 
shaking,  the  mixture  was  centrifuged  at  4300  x  g  for  10  min  at 
5°C  (Allegra  25R,  Beckman  Coulter,  Palo  Alto,  CA),  and  the  su¬ 
pernatant  was  discarded.  A  second  wash  was  performed  with 
the  addition  of  6-7  mL  acetone  to  the  packed  pellet  material, 
following  which  the  samples  were  placed  on  a  shaker  for  10 
min.  The  samples  were  again  centrifuged  (4300  x  g  for  10 
min  at  5°C)  and  the  supernatant  discarded.  This  second  ace¬ 
tone  wash  was  repeated  for  a  total  of  three  washes.  The  pre 
cipitated  protein  samples  were  placed  in  a  desiccator  and 
allowed  to  thoroughly  dry  at  room  temperature.  Following 
drying,  200  pL  of  1M  NaOH  was  added  to  approximately  25  mg 
of  protein.  The  mixture  was  heated  to  70°C  until  the  solid  ma¬ 
terial  dissolved  (approximately  2  h).  Then  20  pL  of  TDG-d8  (300 
ng/mL)  was  added  and  mixed  thoroughly.  Following  the  addi¬ 
tion  of  70  pL  of  3M  HC1, 1.0  mL  ethyl  acetate  was  added  to  the 
samples,  and  then  they  were  vortex  mixed  for  10  min.  The 
samples  were  then  dried  with  the  addition  of  200  mg  sodium 
sulfate,  vortex  mixed  for  10  min,  and  centrifuged.  An  aliquot 
(500  pL)  of  the  ethyl  acetate  supernatant  was  removed,  and  200 
mg  of  sodium  sulfate  was  added  for  further  drying.  To  this 
aliquot  of  ethyl  acetate  were  added  15  pL  of  pyridine,  20  pL  of 
PFB  chloride,  and  10  pL  of  methanol  to  form  the  PFB  deriva¬ 
tives  of  TDG  and  TDG-d8.  The  mixture  was  then  dried  over 
sodium  sulfate  (200  mg).  Bond  Elut  silica  100  mg  Si  SPE  car¬ 
tridges  were  conditioned  with  1  mL  of  ethyl  acetate  before 
use.  A  200-pL  portion  of  the  reacted  mixture  (ethyl  acetate 
containing  the  derivatized  sample)  was  passed  over  the  pre¬ 
conditioned  cartridge.  Finally,  a  400-pL  volume  of  clean  ethyl 
acetate  was  added  to  the  cartridge.  The  two  fractions  (200  and 
400  pL)  were  combined  for  the  GC-MS  analysis. 

Calibration  curves 

Calibration  curves  were  prepared  in  human  and  rat  plasma 
at  six  concentration  levels  of  HD  (12.5,  25,  50, 100,  200,  and 
400nM).  Calibration  curves  prepared  from  human  plasma  were 
made  on  five  separate  days,  whereas  those  from  rat  plasma 
were  made  on  six  different  days.  For  all  samples,  sulfur  mus¬ 
tard  in  saline  (40,000nM,  100  pL)  was  spiked  into  plasma  (9.9 
mL)  to  give  a  final  concentration  of  400nM.  This  sample  was 
serially  diluted  with  blank  plasma  (containing  no  1  ID)  to  obtain 
the  desired  concentration  levels.  The  plasma  spiked  with  HD 
was  then  incubated  for  2  h  at  40°C.  Protein  precipitated  from 
the  calibration  curve  samples  was  processed  as  described  in 
Sample  preparation.  The  actual  HD  concentration  calculated 
for  each  day  was  derived  from  its  respective  calibration  curve. 
Across  days,  the  mean  concentration,  standard  deviation  (SD), 
percent  coefficient  of  variation  [ % CV;  (SD/mean)  x  100]  and 
percent  (%)  error  (percent  difference  between  calculated  and 
expected  concentrations)  were  calculated. 


Intraday  and  interday  precision  and  accuracy  studies 

Intraday  and  interday  variability  were  used  to  assess  preci¬ 
sion  and  accuracy  of  the  method.  For  these  studies,  test  un¬ 
known  samples  were  made  at  three  HD  levels  (25,  150,  and 
300nM)  to  represent  low,  medium,  and  high  concentrations  on 
the  calibration  curve.  To  prepare  the  samples,  human  plasma 
was  spiked  to  a  level  of  400nM  HD.  This  standard  was  serially 
diluted  with  plasma,  resulting  in  HD  standard  concentrations 
of  0,  25,  50, 100,  200,  and  400nM.  Equal  volumes  of  the  400 
and  200nM  standards  were  mixed  to  obtain  the  300nM  test  un¬ 
known;  similarly  equal  volumes  of  the  100  and  200nM  were 
mixed  to  make  the  150nM  test  unknown.  For  intraday  studies, 
five  sets  (n  =  5)  of  test  unknown  samples  were  prepared  and 
analyzed  together.  Precision  was  assessed  by  expressing  the 
%CV.  Accuracy  was  expressed  as  %  error.  Interday  test  un¬ 
known  samples  were  prepared  as  described  except  that  the 
samples  were  prepared  and  analyzed  on  five  different  days.  A 
calibration  curve  was  prepared  and  analyzed  prior  to  each 
(intra-  and  interday)  study.  All  test  unknown  samples  were 
quantified  by  using  the  mean  values  for  these  calibration 
curves.  The  lOOnM  standards  from  the  test  unknown  prepara¬ 
tion  scheme  described  were  utilized  as  check  standards  to 
monitor  assay  performance.  Four  check  standards  were  ana¬ 
lyzed  during  the  intraday  study.  For  the  interday  study,  a  total 
of  nine  check  standards  were  run,  one  at  the  end  of  day  1  and 
thereafter  at  the  beginning  and  end  of  each  subsequent  study 
day  (days  2  5). 

Instrumentation 

GC-MS  separations  were  performed  on  an  Agilent  6890  GC 
interfaced  to  an  Agilent  5973  mass-selective  detector.  The  GC 
was  fitted  with  a  30-m  x  0.25  mm  internal  diameter  DB-5MS 
bonded  phase  capillary  column  (0.25-pm  film  thickness,  J&W 
Scientific).  Helium  was  used  as  the  carrier  gas  in  the  con¬ 
stant  flow  mode  (1.2  mL/min).  The  oven  temperature  was  held 
initially  at  80°C  for  1  min,  programmed  from  80  to  225°C  at 
30°C/min,  and  held  at  225  'C  for  11  min.  Splitless  injections  of 
1  pL  volume  were  made  usrng  an  Agilent  7683  autosampler. 
The  injection  port  temperature  was  set  at  250  C,  split  vent 
delay  at  1  min,  and  the  transfer-line  temperature  at  280°C.  Typ¬ 
ical  retention  times  were  13.3  min  for  TDG-d8PFB  and  13.5 
min  for  TDGPFB. 

The  MS  analysis  was  conducted  using  NICI  with  methane  as 
the  reagent  gas.  The  resolution  capability  of  the  Agilent  5973 
mass  selective  detector  is  specified  as  unit  mass.  The  source 
and  quadrupole  temperatures  were  set  at  150  and  120  C,  re¬ 
spectively.  The  molecular  ion  and  the  M  + 1  isotope  ion  were 
monitored  for  the  TDGPFB  (m/z  510  and  511)  and  TI)G-d8PFB 
internal  standard  (m/z  518  and  519).  Relative  to  monitoring  a 
single  ion,  two  ions  were  used  to  provide  further  confirmation 
of  peak  identity.  The  m/z  511  and  m/z  519  were  utilized  strictly 
as  qualifier  ions.  However,  single  ions  for  both  analyte  TDGPFB 
(m/z  510)  and  internal  standard  TDG-d8PFB  (m/z  518)  were 
used  for  quantification.  The  dwell  time  for  each  ion  was  100 
ms,  resulting  in  a  total  scan  rate  of  2.15  cycles/s.  The  electron 
multiplier  voltage  was  set  at  +400  volts  relative  to  the  autotune 
setting. 
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Data  analysis 

Graphic  representations,  linear  regression  analyses,  calcu¬ 
lations  for  HD  concentrations,  and  other  descriptive  statistics 
were  accomplished  with  GraphPad  Prism  software  (version 
3.02,  GraphPad  Software,  San  Diego,  CA)  or  Microsoft  Excel 
2000.  Ions  used  for  quantification  purposes  were  m/z  518  and 
519  for  the  internal  standard  (TI)G-d8PFB)  and  m/z  510  and 
511  for  the  analyte  (TDGPFB).  The  area  under  the  curve  (AFC) 
was  determined  for  TDGPFB  and  TDG-d8PFB  and  was  ex¬ 
pressed  as  an  AUC  ratio  (TDGPFB/TDG-d8PFB),  which  was 
then  normalized  by  the  mass  of  protein  utilized  for  each  sample 
(AUC  ratio/mg  protein).  Calibration  curves  depicting  the  rela¬ 
tionship  between  AUC  ratio/mg  protein  as  a  function  of  HD 
plasma  concentration  (nM)  were  utilized  to  express  data  as 
HD  plasma  concentration  (nM). 


Results  and  Discussion 

The  assay  reported  in  this  study  has  been  utilized  previously 
in  research  animal  models  and  in  an  accidental  human  expo¬ 
sure  to  confirm  HD  as  the  agent  involved  (13,16).  The  purpose 
of  the  study  reported  herein  is  1.  to  characterize  the  perfor¬ 
mance  of  the  assay  in  terms  of  variability,  linearity,  and  sensi¬ 
tivity  and  2.  to  use  the  method  to  confirm  delivery  of  different 
doses  of  HD  administered  as  an  inhaled  nebulized  aerosol  from 
the  lung  to  the  blood  compartment  in  a  rat  exposure  model. 

Calibration  curves 

Calibration  curves  were  prepared  from  commercially  pur¬ 
chased  human  and  rat  plasma.  The  purpose  was  to  examine  the 
linearity  and  variability  (as  determined  by  %CVand  %  error)  in 


Table  i.  Calibration  Curves  from  Human  Plasma 
Individual  Days* 

on 

Sulfur  Mustard  Standard  (Expected  Concentration! 

12.5 

25 

50 

100 

200 

400 

nM 

nM 

nM 

nM 

nM 

nM 

Day 

Calculated  (nM) 

1 

6.48 

24.4 

50.0 

107.4 

213.3 

391.7 

2 

8.39 

26.6 

51.2 

109.0 

203.5 

395.9 

3 

8.25 

21.2 

51.9 

109.1 

206.2 

394.8 

4 

9.67 

27.2 

50.1 

106.0 

201.8 

397.5 

5 

8.88 

24.2 

51.6 

106.6 

206.7 

394.9 

Mean 

concentration 

8.33 

24.7 

51.0 

107.6 

206.3 

395.0 

(nM) 

SD 

1.17 

2.36 

0.87 

1.39 

4.37 

2.11 

%CV 

14.1 

9.55 

1.71 

1.30 

2.11 

0.  53 

%  Error 

-33.3 

-1.20 

2.00 

7.60 

3.15 

-1.25 

*  The  calculated  HD  concentrations  were  derived  from  respective  daily  calibration 
curves.  The  %  coefficient  of  variation  (%CV)  was  calculated  by  (SD/mean)  x  100 
(SD  =  standard  deviation).  The  %  error  was  calculated  by  percent  difference 
between  calculated  and  expected  concentrations. 

the  matrices  from  both  species.  Linear  regression  analysis  for 
human  plasma  samples  indicated  that  the  range  of  correlation 
coefficients  (R2)  was  0.9974-0.9998  for  individual  days  with  a 
value  for  the  mean  data  of  0.9987.  Similarly  for  rat  plasma,  the 
range  for  individual  days  was  0.9990  to  0.9999  with  the  value 
for  the  mean  data  of  0.9992.  Calibration  curves  constructed 
from  human  and  rat  plasma  data  demonstrated  consistency 
from  day  to  day.  The  %CVs  derived  from  the  mean  concentra¬ 
tion  data  at  levels  12.5,  25.0,  50.0, 100,  200,  and  400nM  were 
14.1, 9.55, 1.71, 1.30, 2.11,  and  0.53%,  respectively,  in  human 
(n  =  5)  and  10.63,  7.55,  2.34,  6.19,  2.65,  and  0.57%,  respec¬ 
tively,  in  rat  ( n  =  6)  plasma  (Tables  I  and  II).  The  data  indicate 
that  the  %  errors  were  in  the  acceptable  range  in  accordance 
with  FDA  guidance  (i.e.,  less  than  15%)  (17),  except  for  the 
12.5nM  level  in  human  plasma,  which  gave  a  mean  concen 
tration  of  8.33nM,  resulting  in  an  error  of  33.3%.  Data  from  the 


Table  II.  Calibration  Curves  from  Rat  Plasma  on 

Individual  Days* 

Sulfur  Mustard  Standard  (Expected  Concentration) 

12.5 

25 

50 

100 

200 

400 

nM 

nM 

nM 

nM 

nM 

nM 

Day 

Calculated  (nM) 

1 

11.44 

23.16 

49.00 

108.26 

199.56 

398.43 

2 

11.33 

26.47 

50.74 

95.05 

209.15 

396.51 

3 

12.15 

26.64 

49.93 

99.91 

201.06 

399.41 

4 

14.43 

29.11 

49.02 

96.28 

196.52 

402.47 

5 

13.82 

28.07 

50.69 

90.04 

208.18 

398.08 

6 

13.90 

26.78 

47.75 

98.14 

197.83 

401.68 

Mean 

concentration 

12.85 

26.71 

49.52 

97.95 

202.05 

399.43 

(nM) 

SD 

1.37 

2.03 

1.16 

6.06 

5.36 

2.26 

%  CV 

10.63 

7.55 

2.34 

6.19 

2.65 

0.  57 

%  f  rror 

2.77 

6.82 

-0.96 

-2.05 

1.03 

-0.14 

*  The  calculated  HD  concentrations  were  derived  from  respective  daily  calibration 
curves.  The  %  coefficient  of  variation  (%CV)  was  calculated  by  (SD/mean)  x  100 
(SD  =  standard  deviation).  The  %  error  was  calculated  by  percent  difference 
between  calculated  and  expected  concentrations. 

Table  III.  Precision  and  Accuracy  for  Intraday  Test 
Samples  in  Human  Plasma* 


Expected 

C  oncentration 
(nM) 

Mean  (n  =  5) 
Calculated 
f  oncentration 
(nM) 

Precision 

Accuracy 
%  Error 

SD 

%CV 

300.0 

292.94 

5.52 

1.88 

-2.35 

150.0 

153.59 

21.18 

13.79 

2.39 

25.0 

22.03 

1.29 

5.86 

-11.87 

0.0 

0 

0 

0 

0 

*  Five  (N  5)  sets  of  test  samples  (0,  25,  1 50,  and  300nM)  in  human  plasma  were 
made  at  once  and  analyzed  together.  Precision  was  assessed  by  expressing  the 
%  coefficient  of  variation  (%CV)  calculated  by  'SD/mean)  x  100  (SD  =  standard 
deviation).  Accuracy  was  expressed  as  %  error  and  was  calculated  by  percent 
difference  between  calculated  and  expected  concentrations. 
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calibration  curves  in  human  and  rat  plasma  demonstrate  that 
the  method  is  reproducible  across  all  concentration  levels  em 
ployed  in  the  study  as  indicated  by  the  %CVs,  which  were  less 
than  15%.  However,  accuracy  at  the  low  level  (12.5nM)  in  the 
human  data  is  relatively  large  (33.3%).  This  is  because  the 
value  is  below  the  limit  of  Quantitation  (LOQ)  (20.88nM)  de 
termination  (see  Limits  of  detection  and  quantitation). 

Precision  and  accuracy 

Studies  to  examine  intra-  and  interday  variability  were  as¬ 
sessed  by  calculating  precision  and  accuracy  values  for  test 
samples  in  human  plasma.  The  concentration  levels  (25, 150, 
300nM)  used  in  the  studies  were  chosen  to  represent  low, 
medium  and  high  concentrations  of  HD  relative  to  those  em¬ 
ployed  in  the  calibration  curve.  Data  demonstrating  precision 
and  accuracy  of  intra-  and  interday  test  samples  are  shown  in 


Table  \X.  Precision  and  Accuracy  for  Interday  Test 
Samples  in  Human  Plasma* 


Expected 

Concentration 

(nM) 

Mean  (n  =  5) 
Calculated 
Concentration 
(nM) 

Precision 

Accuracy 
%  Error 

SD 

%CV 

300.0 

307.60 

12.14 

3.95 

2.53 

150.0 

153.93 

4.10 

2.66 

2.62 

25.0 

26.18 

3.83 

14.64 

4.73 

0.0 

0 

0 

0.00 

0.00 

*  Five  (N  =  5)  sets  of  test  samples  (0,  25,  ]  50  and  300nM)  in  human  plasma  were 
made  and  analyzed  daily  for  five  days.  Precision  was  assessed  by  expressing  the 
%  coefficient  of  variation  (%CV)  calculated  by  (SD/mean)  x  1 00  (5D  =  standard 
deviation).  Accuracy  was  expressed  as  %  error  and  was  calculated  by  percent 
difference  between  calculated  and  expected  concentrations. 


Table  V.  Check  Standards* 

Intraday 

Interday 

Standard 

Calculated 

Study 

Calculated 

number 

concentration  (nM) 

day 

concentration  (nM) 

1 

87.67 

1 

94.81 

2 

96.72 

2 

97.88 

3 

89.07 

2 

93.30 

4 

9S.67 

3 

103.24 

- 

- 

3 

105.94 

- 

- 

4 

100.51 

- 

- 

4 

102.92 

- 

- 

5 

94.23 

- 

- 

5 

97.47 

Mean 

92.28 

98.92 

SD 

4.57 

4.47 

%CV 

4.96 

4.52 

%  Error 

-7.72 

-4.18 

*  Check  standards  made  at  a  level  of  1  OOnM  were  run  to  monitor  assay 
performance  during  the  course  of  intra-  and  interday  studies.  For  the  intraday 
study,  four  check  standards  were  run  during  the  course  of  the  study  day.  For 
interday  studies,  the  check  standards  were  run  at  the  end  of  day  1  and  thereafter 
at  the  beginning  and  end  of  each  subsequent  study  day  (interday  2-5)  for  a  total 
of  nine  check  standards. 

Tables  III  and  IV,  respectively.  For  intraday  samples,  %CVs 
were  5.86,  13.79,  and  1.88  for  concentrations  25,  150,  and 
300nM,  respectively.  The  accuracy  as  reflected  by  %  error  was 
-11.87, 2.39,  and  -2.35%  for  25, 150,  and  300nM,  respectively. 
The  %CVs  for  the  interday  samples  were  14.64, 2.66,  and  3.95 
and  %  errors  were  4.73,  2.62,  and  2.53%  for  25,  150,  and 
300nM,  respectively.  These  data  indicate  that  the  assay  demon¬ 
strates  precision  and  accuracy  at  the  studied  HD  levels,  which 
are  acceptable  as  defined  by  FDA  guidance  (i.e.,  less  than  15%) 
(17).  To  monitor  the  assay  performance,  check  standards 
(lOOnM)  were  run  during  the  course  of  the  intra-  and  interday 
studies.  The  calculated  concentrations  for  the  check  standards 
run  during  the  intra-  and  interday  studies  are  presented  in 
Table  X.  The  mean  calculated  concentration  for  check  standards 
run  during  the  intraday  study  was  92.28nM  with  a  %CV  of 
4.96%  and  %  error  of  -7.72%.  The  calculated  concentrations 
for  the  check  standards  (lOOnM)  run  during  the  interday  study 
are  presented  in  Table  V.  The  mean  calculated  concentration 
was  98.92nM  with  a  %CV  and  %  error  of  4.52  and  -4.18%,  re¬ 
spectively. 

Limits  of  detection  (LOD)  and  quantitation  (LOQ) 

The  LOD  and  LOQ  were  determined  by  the  relationship  3S0 
and  10S0  respectively  (18,19).  The  S()  was  calculated  as  pre¬ 
sented  (19)  and  derived  from  the  //-intercept  of  plots  of  stan¬ 
dard  deviations  from  calibration  curves  as  a  function  of  the 
respective  HD  concentrations  (12.5, 25, 50, 100,  and  200nM). 
For  these  assessments,  data  from  seven  calibration  curves  run 
over  a  five-week  period  were  utilized  from  human  plasma.  For 
rat  studies,  data  from  calibration  curves  run  along  with  un¬ 
known  samples  from  all  the  inhalation  studies  were  utilized. 
The  rat  studies  in  total  encompassed  a  period  of  18  months. 
The  analysis  indicated  that  the  LOD  and  LOQ  for  rat  and 
human  plasma  were  similar.  The  LOD  was  determined  as 
6.26nM  for  human  and  5.02nM  for  rat,  and  LOQ  was  20.88nM 
for  human  and  16.73nM  for  rat  plasma.  The  LOQ  determina¬ 
tion  in  human  plasma  is  consistent  with  the  relatively  large  % 
error  obtained  from  the  lowest  level  (12.5nM,  33.3%)  in  the 
calibration  curves  from  human  plasma.  That  is,  at  that  level  the 
assay  did  not  produce  quantifiably  valid  results.  For  a  similar 
assay,  Capacio  et  al.  (13)  reported  the  lowest  level  on  the  cali¬ 
bration  curve  used  to  quantify  samples  as  25nM,  because  that 
was  the  lowest  standard  utilized  on  the  calibration  curve  with 
a  signal/noise  ratio  of  a  3:1.  Although  the  LOQ  determined  in 
this  study  (20.88nM)  was  derived  differently,  it  is  similar  to  that 
reported  by  Capacio  et  al.  (13).  Recently,  numerous  modifica¬ 
tions  to  the  procedure  described  in  this  present  study  have 
been  implemented  (20).  They  have  resulted  in  approximately  a 
10-fold  increase  in  sensitivity  (LOQ  =  1.6nM)  relative  to  that  re¬ 
ported  in  this  present  study  (20). 

Rat  inhalation  studies 

The  utility  of  the  method  was  demonstrated  in  a  rat  inhala¬ 
tion  model  using  nebulized  HD.  The  graphic  representation  of 
HD  plasma  concentration  as  a  function  of  time  for  each  dose  is 
depicted  in  Figure  1.  In  this  study,  HD  was  detectable  in  the 
plasma  of  the  exposed  animals  out  to  the  last  time  point  ex 
amined  (24  h).  The  data  indicate  a  dose  dependent  relationship 
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between  maximal  HD  plasma  concentrations  and  dose  admin¬ 
istered.  Mean  maximal  plasma  HD  concentrations  were  33.61, 
107.1,  335.6,  390.7,  566.0,  and  513.5nM  for  doses  250,  750, 
1750, 2250, 2750,  and  3000  pg,  respectively.  Regression  anal- 


Figure  1.  Plasma  concentration-time  profile  of  inhaled  nebulized  HD  in 
rats.  Animals  received  HD  directly  into  the  lung  as  a  nebulized  aerosol 
at  six  doses  over  10  min:  2250, 2750,  and  3000  pg  HD  (A)  and  250,  750, 
and  1750  pg  HD  IB).  Blood  was  obtained  from  4-6  animals  per  time 
point  at  0,  0.5  1,  3,  fa,  and  24  h  after  exposure.  Data  are  mean  HD 
plasma  concentration  plus  or  minus  the  standard  deviation. 
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Figure  2.  Maximal  plasma  concentration  of  HD  and  dose  administered. 
Animals  nebulized  HD  aerosol  for  total  doses  of  250,  750, 1 750, 2250, 
2750,  and  3000  pg.  The  maximal  HD  plasma  concentration  obtained  (± 
SDi  as  a  function  of  dose  delivered  are  plotted.  The  data  indicate  a 
linear  relationship  as  reflected  by  a  correlation  coeffic  ient  of  0.9726. 


ysis  of  mean  maximal  plasma  HD  concentration  as  a  function 
of  dose  produced  an  R 2  of  0.9728  (Figure  2).  Although  the  re¬ 
gression  analysis  indicates  a  linear  relationship  between  dose 
and  maximal  plasma  concentration,  the  relationship  does  not 
appear  to  hold  at  the  two  highest  doses  with  respect  to  the 
maximal  plasma  concentrations  and  the  overall  time-course. 
Potential  explanations  for  the  observation  include  animal  to 
animal  variability  and/or  saturation  of  processes  in  the  lung 
that  are  involved  in  partitioning  of  HD  into  the  circulation  at 
these  higher  levels  of  exposure.  More  studies  are  necessary  to 
further  elucidate  the  dose-plasma  concentration  relationship  at 
this  level  of  HD  administration. 

This  is  the  first  study  to  examine  the  systemic  presence  of  an 
HD  biomarker  from  inhaled  and  nebulized  HD  driven  directly 
into  the  lung  by  ventilatory  action  in  the  anesthetized  rat.  In 
other  HD  inhalation  models  (i.e.,  nose  only  administration) 
(21),  many  of  the  systemic  effects  may  not  be  observed  since 
detoxification  processes  take  place  in  the  nasal  passages.  This 
phenomenon  can  scrub  a  significant  amount  of  an  inhaled 
toxicant  in  obligate  nose  breathers  (22),  and  data  regarding  the 
true  effects  of  inhaled  HD  on  lung  toxicity  and  the  role  of  the 
lung  in  systemic  absorption  processes  may  not  be  clear.  In  ad¬ 
dition,  the  well-known  effect  of  airway  mucus  plug  formation 
can  lead  to  ingestion  and  subsequent  systemic  HD  poisoning 
via  absorption  from  the  gut.  Models  using  whole-body  vapor  ex¬ 
posures  of  other  chemical  agents  do  not  characterize  systemic 
absorption  exclusively  through  the  lung,  but  represent  a  com 
posite  of  numerous  processes  that  contribute  to  total  blood 
levels  (23).  As  such,  these  models  do  not  truly  represent  HD 
that  enters  the  systemic  circulation  entirely  via  the  lung.  Ad¬ 
vantages  to  the  model  utilized  in  the  experiments  reported 
herein  are  that  it  enables  the  study  of  HD  that  is  directly  fo¬ 
cused  on  the  lung  and  it  minimizes  interactions  with  other  sys¬ 
tems  such  as  airway  components  or  skin.  Therefore,  a  direct 
relationship  between  dose  delivered  and  blood  levels  and  or  tox¬ 
icity  can  be  obtained.  The  model  will  provide  significant  data 
regarding  the  true  effects  of  inhaled  HD  and  the  role  of  the 
lung  in  systemic  absorption,  distribution,  metabolism,  and  ex¬ 
cretion  processes. 


Conclusions 

A  previously  reported  GC-MS  method  for  determining  in 
vivo  exposure  to  III)  has  been  validated.  The  method  leverages 
the  binding  of  HD  to  plasma  proteins  with  subsequent  cleavage 
and  derivatization  of  the  adduct  and  GC-MS-NICI  analysis. 
Results  from  this  study  indicate  an  accurate,  precise  and 
sensitive  method.  Variability  determined  in  precision  and 
accuracy  studies  was  less  than  15%.  The  LOQs  determined  in 
human  and  rat  plasma  were  calculated  to  be  20.88nM  and 
16.73nM,  respectively.  The  use  of  the  method  in  a  rat  inhala¬ 
tion  model  demonstrates  the  presence  of  the  plasma  protein 
biomarker  in  a  dose-dependent  fashion  after  the  delivery  of 
nebulized  HD  to  the  lung.  Regression  analysis  of  mean  max¬ 
ima1  plasma  HD  concentration  as  a  function  of  dose  produced 
an  R2  of  0.9728. 
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